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Monopole Condensation in Lattice SU(2) QCD * 
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Department of Physics, Kanazawa University, Kanazawa 920-11, Japan 

ABSTRACT 

This is the short review of Monte-Carlo studies of quark confinement in lattice 
QCD. After abelian projections both in the maximally abelian and Polyakov 
Q\ \ gauges, it is seen that the monopole part alone is responsible for confinement. A 

block spin transformation on the dual lattice suggests that lattice SU (2) QCD 
is always ( for all (3) in the monopole condensed phase and so in the confinement 
^ \ phase in the infinite volume limit. 
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It is crucial to understand the mechanism of quark confinement in order to explain 
hadron physics out of QCD. Our standpoint is based on the 'tHooft idea of abelian 
projection of QCD.0 The abelian projection is to fix the gauge in such a way that 
the maximal torus group remains unbroken. After the abelian projection, monopoles 
appear as a topological quantity in the residual abelian channel. QCD is reduced to 
an abelian theory with electric charges and monopoles. If the monopoles make Bose 
condensation, charged quarks and gluons are confined due to the dual Meissner effect. 

Based on this standpoint, the present author and his collaborators have studied 
color confinement mechanism and hadron physics performing Monte Carlo simulations 
of abelian projection in lattice QCD.ml3 The aim of the study is to ascetain correctness 
^ ■ of the picture, that is, to check if monopole condensation really occurs in QCD. Here 

I review the results of these studies compactly.^ 



2. Abelian dominance and monopole dominance 

2.1. Maximally abelian and Polyakov gauges 

There are infinite ways of abelian projection extracting such an abelian theory 
out of QCD. We have found two gauges which show interesting behaviors called abelian 
dominanceSi'0 



*This is based on works done in collaboration with Y.Matsubara, H.Shiba, S.Kitahara, S.Ejiri, S.Okude 
and K.Yotsuji. 

^E-mail address: suzuki@hep.s. kanazawa-u.ac.jp 

* We have also studied hadron physics based on andaftared effective Lagrangian constructed directly 
from QCD on the assumption of the above picture,Ejc2l but this time I skip all of them. 
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SU(2) POLYAKOV GAUGE 
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Fig.l Monopole Dirac string 
and photon contributions to 
Polyakov loops in MA gauge. 



Fig. 2 Monopole Dirac string and 
photon contributions to Polyakov 
loops in the Polyakov gauge. 



One is the maximally abelian (MA) gauge. Define a matrix in 577(2) QCD 

x ( s ) = I][^(' s '/ u ) or 3^7 t (s, / u) + U j (s - fafa)a 3 U(s - p,,/i)] 

= Xi(s)<7i + X 2 (s)a 2 + X 3 (s)a 3 . 

Then a gauge satisfying Xi(s) = X 2 (s) = is the MA gaug which tends to a U(l)- 
covariant gauge (d^ ^zigA'^ L )A ± ^ = in the continuum limit. The other is the Polyakov 
gauge which is defined by diagonalizing Polyakov loop operators. We have found that 
important features of confinement, i.e., the string tension and the charateristic be- 
haviors of the Polyakov loops are well reproduced in terms of the abelian operators 
0(u(s, fa)) both in MA gaugeii and in the Polyakov gauge.0 See Figures 1 ~ 4. 

2.2. Monopole dominance in the Polyakov loops 

An abelian Polyakov loop P which is written in terms of abelian link fields alone 
is given by a product of contributions from Dirac strings of monopoles and from pho- 
tons.00 We have observed the photon and the Dirac-string contributions separately. 
See Figs.l and 2. The characteristic features of the Polyakov loops as an order param- 
eter of the deconfinement transition are due to the Dirac string contributions alone. 
The photon part has a finite contribution for both phases and it changes only slightly. 
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Fig. 3 Monopole and photon 
contributions to the string 
tension in MA gauge in SU(2) 

QCD in comparison with the full Fi S- 4 Polyakov loop correlations 

and the abelian ones. in the Polyakov gauge. 



The fact that monopoles are responsible for the essential feature of the Polyakov loop 
is found also in U(l) and in SU(3). It is remarkable that the behavior of the abelian 
Polyakov loops as an order parameter and the monopole responsibility are seen in any 
gauged This is the first phenomena suggesting gauge independence of the 'tHooft 
conjecture. 

2.3. The string tension and monopoles 

Shiba and SuzukiBi have shown furthermore that monopoles alone can reproduce 
the full value of the string tension in SU(2) QCD. An abelian Wilson loop operator 
is rewritten by a product of monopole and photon contributions. The string tensions 
evaluated from monopole and photon contributions are plotted in Fig. 3. The string 
tension is reproduced only from the monopole contribution. 

It is possible to know the static potential from the correlation of the Polyakov 
loops. The string tension is then derived from the static potential. We have measured 
the correlations using the Polyakov loop operators written in terms of full and abelian 
link fields and link fields given by the monopole Dirac string. The preliminary data are 
plotted in Fig. 4 in the case of the Polyakov gauge. A linear behavior is seen with the 
almost similar string tension in all three cases. This is another data supporting gauge 
invariance of confinement mechanism due to monopole condensation. 
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3. Monopole action and condensation 

The above abelian dominance suggests that a set of U(l) invariant operators 
{0(u(s, /j,))} are enough to describe confinement. Then there must exists an effective 
U{\) action S e jf(u) describing confinement. We tried to derive S e ff(u) using Schwinger- 
Dyson equations, but failed to get it in a compact and local form.B S e ff(u) contains 
larger and larger loops as (3. 

Shiba and Suzuki author tried to perform a dual transformation of S e ff(u) in 
SU(2)QCB and to obtain the effective U(l) actioniffl in terms of monopole cur- 
rents.^ To study the lone range behavior important in QCD, They have considered 
also extended monopoles.o The extended monopole currents are defined by the number 
of the Dirac strings surrounding an extended cube: 

n-l 

k ( ^(s) = E k^ns + {n-l)p, + iC>+j P + la), (1) 

i,j,l=0 

where fc^(s) is the ordinary monopole current.0 Considering extended monopoles cor- 
responds to performing a block spin transformation on a dual latticeUHl and so It is 
suitable for exploring the long range property of QCD. 

The partition function of interacting monopole currents is expressed as 

oo 

Z =(R E )(U S ^ is ) fi )eM-S[k)). (2) 

s,/t fc M (s)=-oo s 

It is natural to assume S[k] = ^2ifiSi[k\. Here /j is a coupling constant of an in- 
teraction Si[k}. For example, fi is the coupling of the self energy term J^n^ik^s)) 2 , 
f 2 is the coupling of a nearest-neighbor interaction term J2 n ,/i k^(s)k^(s + fi) and f 3 
is the coupling of another nearest-neighbor term J2 n ,^k ll (s)k f j - (s + z>).0 Shiba and 
SuzukiBi'l extended a method developed by SwendsenM to the system of monopole 
currents obeying the current conservation rule. 

The monopole actions are obtained locally enough for all extended monopoles con- 
sidered even in the scaling region. They are lattice volume independent. The coupling 
constant /i of the self-energy term is dominant and the coupling constants decrease 
rapidly as the distance between the two monopole currents increases. 

To study monopole dynamics, we have also studied the length of monopole loops 
and monopole charge distribution. We have found that l)each vacuum has only one 
long connected loop and a small number of short loops. The difference of the two types 
of loops is very clear. 2)The length of the long loop becomes shorter as (3 becomes 
larger. In the deep deconfinement region, there disappears a long loop. 3)Monopole 
charges are almost ±1. 

Since the long loop is important, we plot the value of monopole action versus 
length of the long loops as seen in Fig. 5. The value of the action is proportional to the 
length L of the long loop is well approximated by j\ x L. 
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Action versus monopole length 
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Fig. 5 Monopole action versus 
length of loops 
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Fig. 6 Coupling constants fi 
versus (3 for 2 3 , 3 3 , and 4 3 
monopoles on 24 4 lattice. 



As done in compact QED,c3 the entropy of a monopole loop can be estimated 
as In 7 per unit loop length. Since the action is approximated by the self energy part 
fiL, the free energy per unit monopole loop length is approximated by (j\ — In 7) 
. If fx < In 7, the entropy dominates over the energy, which means condensation of 
monopoles. In Fig. 6, f\ versus (3 for various extended monopoles on 24 4 lattice is 
shown in comparison with the entropy value In 7. Each extended monopole has its own 
(3 region where the condition fi < In 7 is satisfied. When the extendedness is bigger, 
larger (3 is included in such a region. Larger extended monopoles are more important 
in determining the phase transition point. 

The behaviors of fi are different for different extended monopoles. However, if 
we plot them versus b = n x a((3), we get a unique curve as in Fig. 7. The coupling 
constants seem to depend only on b, not on the extendedness nor (3. There is a critical 
b c corresponding to critical i.e., b c = na(/3"). 

Now we can derive important conclusions. Suppose the effective monopole action 
remains the same for any extended monopoles in the infinite volume limit. Then the 
finiteness of b c = na((3™) suggests /3" becomes infinite when the extendedness n goes 
to infinity. SU(2) lattice QCD is always (for all j3) in the monopole condensed and 
then in the color confinement phaseS This is one of what one wants to prove in the 
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Fig. 7 Coupling constants fi 
and fi versus b. 
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Fig. 8 Renormalization flow in 
the fi — fi plane 



framework of lattice QCD. 

Notice again that considering extended monopoles corresponds to performing a 
block spin transformation on the dual lattice. The above fact that the effective actions 
for all extended monopoles considered are the same for fixed b means that the action 
may be the renormalized trajectory on which one can take the continuum limit. See 
Fig.8. Our results suggest the continuum monopole action takes the form predicted by 
Smit and Sijs.@ 

4. Simulation of monopole action 

To test if the fixed monopole action is the renormalized trajectory, we have tried 
to simulate the system with the action.El 

We have obtained the following results: 

1. The action with quadratic interactions alone is not good, since it gives a longer 
monopole loop and a larger string tension than expected. We need a term which 
gives a repulsive force between monopole currents. 

2. With such small repulsive terms, we get the nice fit of the loop length and the 
string tension on small lattices. Finite-size effects are seen to be small. 

3. We have also made Monte-Carlo simulations in the T ^ QCD. Only the self- 
coupling term can reproduce the qualitative features of the deconfinement tran- 
sition as seen in Fig. 10. 
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the lattice size 



Fig. 10 Polyakov loop from the 
self-coupling monopole 
interactions alone in the 
two-current case 



5. Discussions and outlook 

Some comments are in order. 

1. It is important to find an order parameter of confinement due to monopole con- 
densation.0 The number of Dirac string nn(s) may be a candidate!!! n i4: (s) is a 
magnetic U(l) variant quantity.0 

2. Gauge independence should be proved if the monopole condensation is the real 
confinement mechanism. The data of Polyakov loops from monopoles are encour- 
aging. That the correct string tension could be derived both in MA gauge and 
Polyakov gauge is also encouraging. Gauge independent results will be obtained 
if we go to large (3 on larger lattice. 

3. How to test the correctness of this idea? The theory predicts the existence of 
an axial vector glueball-like state C(J pc = V ) and a scalar glueball-like state 
x(J pc — ++ ). The masses seem to satisfy m c ~ m x .&0 The masses could not be 
too heavy. They have to exist under 2Gev. To evaluate the correlation between 
the state and the light hadrons in MC simulations of full QCD is very important 
to derive the total width and the branching ratios. 
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